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Abstract —In this paper, we analyze dynamic switching net¬ 
works, wherein the networks switch arhitrarily among a set 
of topologies. For this class of dynamic networks, we derive 
an epidemic threshold, considering the SIS epidemic model. 
First, an epidemic prohahilistlc model is developed assuming 
independence between states of nodes. We Identify the conditions 
under which the epidemic dies out hy linearizing the underlying 
dynamical system and analyzing its asymptotic stability around 
the origin. The concept of joint spectral radius is then used 
to derive the epidemic threshold, which is later validated using 
several networks (Watts-Strogatz, Barabasi-Albert, MIT reality 
mining graphs. Regular, and Gilbert). A simplified version of 
the epidemic threshold is proposed for undirected networks. 
Moreover, in the case of static networks, the derived epidemic 
threshold is shown to match conventional analytical results. Then, 
analytical results for the epidemic threshold of dynamic networks 
are proved to be applicable to periodic networks. For dynamic 
regular networks, we demonstrate that the epidemic threshold is 
identical to the epidemic threshold for static regular networks. 
An upper bound for the epidemic spread probability in dynamic 
Gilbert networks is also derived and verified using simulation. 

Index Terms —Dynamic Networks, Epidemic Threshold, Dy¬ 
namical System. 


I. Introduction 

Epidemics typically start with some initial infected nodes. 
There is a probability that healthy neighboring nodes, close to 
the infected nodes, become infected. With time and, in some 
cases, with external intervention, infected nodes recover and 
revert to a healthy state. The study of epidemic dispersal on 
networks aims at understanding how epidemics evolve and 
spread in networks. When an infection enters a network, it 
is very useful to be able to determine whether it will die 
out or become a massive out break. The epidemic threshold 
addresses this question, taking into account both the network 
topology and the epidemic strength. The spread of epidemics 
in static networks has been studied extensively m-iiii. More 
recently, epidemic dispersal in dynamic networks has garnered 
attention |9l- llMI . In a dynamic network, links between nodes 
are functions of time. This dynamic nature presents a more 
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realistic picture of the spread of epidemic infection. Let us 
examine a few examples of dynamic networks which can be 
modeled. First in the area of human epidemics, the nodes 
(human beings) are constantly moving from one location to 
another, thus the contact networks between people change 
over time. In the case of a mobile ad hoc network, each 
device can move independently and change its links to other 
devices frequently, which results in dynamic topologies. Many 
Bluetooth devices are becoming susceptible to viruses such as 
Cabir or Comm Warrior. Another classic example of dynamic 
networks is the networks for the spread of diseases among 
animals or plants where the factors that influence the spread 
of disease-causing spores are typically dynamic ED-Ea. 

In this paper, we consider the Susceptible-Infected- 
Susceptible (SIS) model for epidemic spread. In this model, 
healthy nodes can become infected through infected neighbors; 
infected nodes have a probability of recovery. In the SIS 
model, an infected node after recovery becomes susceptible 
to infection again. We assume that infected nodes have the 
recovery probability <5, and that the probability of an infected 
node infecting its healthy neighbor is /3. There are numerous 
examples in the real world for which the SIS model is the 
right choice to model epidemics. For example, several diseases 
caused by bacteria do not produce immune response in the 
body, thereby allowing the recovered individuals to return to 
the susceptible population. The SIS model can also be used 
for opinion analysis. In a contact network of people, every 
person can have either a positive or negative attitude toward 
a subject and can constantly change from positive to negative 
and vice versa. 

In this paper, we study dynamic switching networks in 
which adjacency matrices are randomly chosen from sets 
of matrices at each step. We do not consider any temporal 
correlation between subsequent adjacency matrices. In other 
words, we assume that the process of choosing the adjacency 
matrix at time index t is independent from adjacency matrices 
chosen at previous time indices. However, this assumption 
does not hold for dynamic networks with strong temporal 
correlations between successive adjacency matrices. 

First, the nonlinear dynamic nature of nodal infection prob¬ 
abilities, assuming independence among the states of nodes, 
will be developed. Then, we will prove that the origin is 
always an equilibrium point of this time-varying dynamical 
system, and its stability depends on network topology and 
values of S and /3. Next, the linearized version of the non¬ 
linear epidemic system is derived to determine asymptotic 
stability of the origin. We show that if the origin is not a 
stable equilibrium of the system, the epidemic spreads. The 


2 


joint spectral radius of a set of matrices is also defined. In 
Theorem [T] we employ the concept of joint spectral radius in 
order to derive the analytical epidemic threshold for dynamic 
networks. In Theorem|2] the simplified version of the epidemic 
threshold for undirected networks is derived. Because the 
epidemic threshold for undirected networks depends only 
on the largest spectral radius of a set of system matrices, 
evaluation of the epidemic threshold is computationally less 
expensive compared to directed networks. In Corollary [T] the 
derived epidemic threshold confirms conventional analytical 
results for static networks, and then the proposed epidemic 
threshold for dynamic networks is extended to periodic net¬ 
works. We also study epidemic spread in dynamic regular 
networks and show that the epidemic threshold for dynamic 
regular networks is identical to that for static regular networks. 
An upper bound for the probability of an epidemic spreading 
in dynamic Gilbert networks is derived. Finally, we simulate 
epidemics in Watts-Strogatz, Barabasi-Albert, Regular, and 
dynamic Gilbert networks in order to validate our analytical 
results. Additionally, we examine our theoretical results in the 
context of real networks by considering MIT reality mining 
graphs 129) . 

The rest of the paper is organized as follows. In Section 
nil we review related prior work on epidemic thresholds in 
dynamic networks. Section [III] contains general analytical 
results for epidemic thresholds for dynamic networks as well 
as simplified epidemic thresholds for special cases of dynamic 
networks. In Section |IVl we use simulation results to validate 
our theoretical analysis. 

II. Related Previous Work 

In mil, the epidemic threshold for an Susceptible-Infected- 
Recovered (SIR) model is derived for a simple class of 
dynamic random networks. In these networks, the number of 
neighbors of a given node is fixed, but its neighbors change 
stochastically as a Poisson process through instantaneous 
neighbor exchanges. Pairs of edges are chosen continually and 
randomly with equal probability, and they are instantaneously 
interchanged. In ma, the authors present a model that de¬ 
scribes an SIS epidemic on dynamic networks using a set of 
ordinary differential equations. The SIS effective degree model 
for a static contact network in Q is modified by introducing 
link activation and deletion rates. The epidemic threshold is 
calculated for this model and shows that the limiting of the 
maximum nodal degree of a network can prevent the outbreak 
of epidemic. In iflSl . the authors study epidemic dynamics on 
an adaptive network, in which susceptible nodes try to avoid 
infection from infected nodes. To achieve this, the susceptible 
nodes cut their links with infected nodes, using a constant 
rewiring probability, and replace them by links with other 
susceptible nodes. The adaptive rewiring increases isolation 
of infected individuals, and simultaneously, contributes to the 
formation of a highly connected susceptible cluster. Conse¬ 
quently, the local effect of rewiring increases the epidemic 
threshold, while the topological effect renders the network 
vulnerable to epidemics. Moreover, the adaptive nature of 
the system leads to the emergence of bistability and limit 


cycles in its dynamical behavior; however, only one continuous 
dynamical transition exists in the static networks. In 13, the 
authors derive the epidemic threshold for dynamic networks 
with alternating (periodic) adjacency matrices. They consider 
the SIS model for epidemic propagation in networks and 
show that if the dynamic behavior of a time-varying network 
can be characterized by T repeating alternating graphs, and 
L = {Ai, A 2 , ■■■, At}, then the system matrix, S, of this 
dynamical system can be expressed as 

T 

i=l 

where the dimension of Aj is n x n (n is the number of 
nodes), I is an n by n identity matrix, and S and /3 denote, 
respectively, the recovery probability and infection probability. 
The authors of 0 prove that if the spectral radius of the 
system matrix is less than 1, the origin is an asymptotically 
stable equilibrium point of the system, and the epidemic 
dies out. This result holds only for cases having repeating 
patterns of adjacency matrices, with the order of repetition 
is preserved. In ifTOl . the authors study malware propagation 
on mobile ad hoc networks. They extend their results for the 
epidemic threshold of periodic networks in 0 to general cases 
in which the repeating order of adjacency matrices can be 
arbitrary. In Theorem I (Mobility model threshold) of their 
paper, they state that if a mobility model can be represented 
as a sequence of connectivity graphs L — {Ai, A 2 , ■■■, At}, 
with one adjacency matrix At for each index t G {1,2,..,T}, 
then the epidemic threshold is 

T = As, (2) 

where As is the largest eigenvalue of the matrix S defined in 
Eq. O. 

This theorem claims that, for a given dynamic network with 
an adjacency matrix arbitrarily chosen from a set of matrices 
at each index, the condition for asymptotic stability is that the 
spectral radius of the matrix S is less than 1. This is different 
from our analytical results in Theorem [T| discussed in Section 
Hm given the same assumptions. 

III. ANALYTICAL RESULTS 

In this section, we develop a dynamical system for epidemic 
spread, assuming spatial independence between states of nodes 
in a given network. A linearized version of the dynamical 
system is then derived to determine the epidemic threshold. 
Next, by employing the joint spectral radius, we quantify the 
epidemic threshold for dynamic networks, proving that the 
epidemic threshold in undirected networks depends only on 
the maximum spectral radius of the set of system matrices. 
Then we extend the results so obtained to static and periodic 
networks. We also show that the epidemic threshold for 
dynamic regular networks is equal to the epidemic threshold of 
static regular networks. Linally, we calculate the upper bound 
for the probability of epidemic spreading in dynamic Gilbert 
networks. 
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Using the above cited assumption, we can write the infection 
probability of each node in the network as 

p,{t + 1) = 1 -p^{t)d - {1 - [1-pj (<)/?], (3) 

where Ni{t) denotes the set of neighbors to node i at index 
t, which is a function of time. Infection probabilities of nodes 
can be interpreted as state variables of a dynamical system. 
Eq. (O shows that infection probabilities at a given index are 
nonlinear functions of infection probabilities of the previous 
index. Therefore, the epidemic dynamical system is nonlinear. 
The corresponding state space of this nonlinear system is the 
subspace [0,1]" in i?", where n is the number of nodes in 
the network. For instance, when n = 2, the state space is a 
rectangle with vertices represented by points of (0, 0), (0,1), 
(1,0), and (1,1) in R^. Given initial infection probabilities 
of the nodes, we can calculate the evolving trajectory of 
infection probabilities in state space. The family of evolving 
trajectories of states in the state space is called the phase 
portrait. The study of the steady state behavior of dynamical 
systems requires hnding the equilibrium points. If P* is an 
equilibrium point, P*(t + 1) = P*{t) = P*. Therefore, we 
can write 

p* = l-p*5-{l-p*) [] [l-p*/3], (4) 

j&Niit) 

where p* is the infection probability of the node in steady 
state, if an asymptotic equilibrium point is present. Eq. @ 
is the equilibrium equation corresponding to node i. In order 
to find equilibrium points of a given epidemic system with 
n nodes, we must solve a system of n equations with n 
unknowns. In the case of dynamic networks, this system 
of equations has been changing with time, and equilibrium 
points, by definition, are static points that satisfy this system 
of equations for all time. An epidemic dynamical system may 
have more than one equilibrium point. From Eq. (|4]i, the origin 
is always an equilibrium point, meaning that, for all values of 
P and 6 and for any arbitrary topologies of network, the origin 
is always an equilibrium point. However, stability of the origin 
depends on /3, 6, and underlying topology of the network. If 
the origin is an asymptotically stable equilibrium point, the 
epidemic dies out. However, the epidemic spreads when the 
origin is an unstable equilibrium point. 

A. Linearization of System Equations 

One way to identify the stability status of an equilib¬ 
rium point of a nonlinear system is to study stability of 
the linearized system at that equilibrium point. In the case 
of epidemic networks, we are interested in determining the 
stability status of the origin. If the origin is an asymptotically 
stable point, the epidemic dies out only if no other equilibrium 
points exist in the subspace [0,1]"; otherwise, asymptotic 
stability is only local. Therefore, we linearize the epidemic 
nonlinear system at the origin. Ignoring the nonlinear terms in 
Eq. 0, we can write 



Fig. 1. Phase portrait of an epidemic network with three nodes in which the 
epidemic dies out. 



Fig. 2. Phase portrait of an epidemic network with three nodes in which the 
epidemic spreads. 

pPt + l) = pPt){l - 5) + ^ Pjit)P. (5) 

We can rewrite 0 in the form of a matrix equation as 

Pt+i = [(1 - S)I + pAt] Pt, (6) 

where Pt = [pi(t)p 2 {t)...pn{t)Y the system state. At is the 
adjacency matrix at index t, and I denotes the n x n identity 
matrix. We define Mt as Mt = [(1 — 5)1 + PA^ the system 
matrix at index t. 

B. Epidemic Threshold 

As mentioned earlier, stability of the origin in the epidemic 
dynamical system determines under what conditions epidemics 
die out. If the origin is an asymptotically stable equilibrium 
point, the system state reaches the origin and infection proba¬ 
bilities of all nodes become zero and remain zero. Fig. [U de¬ 
picts trajectories of state evolution of an epidemic network for 
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different initial infection probabilities. The network has three 
nodes and its adjacency matrix is static. Node 1 is connected to 
nodes 2 and 3, but no link exists between nodes 2 and 3. Each 
axis of the 3D plot in Fig. [T]represents the infection probability 
of one of the three nodes over time. Dots with Pq label in 
state space represent initial infection probabilities, and the 
dot with Peg represents the equilibrium point of the system. 
The probability constants used are 6 = 0.2 and /3 = 0.1. 
We see that all the trajectories reach the origin regardless 
of their initial states. Fig. [T] shows that, for these probability 
constants, the origin is asymptotically stable and epidemics 
die out. In the case of spreading epidemics, the origin is not 
a stable equilibrium point and the state variables converge 
on a non-zero equilibrium point and remain at that point. 
This equilibrium point determines the final fraction of infected 
nodes. Fig. |2] shows trajectories of state evolution of the same 
network depicted in Fig. [T] The initial states are identical to 
initial states in Fig. [T| However, the value of /3 is different. 
In Fig. 121/3 = 0.6, and the epidemic spreads out. For all 
different infection probability vectors, the system state reaches 
the equilibrium point P^q = [0.76791,0.69731,0.69731] with 
elements that represent, respectively, infection probability of 
nodes 1, 2, and 3 in steady state. 

Before tackling the problem of stability of the origin, 
essentials definitions must be presented: 


Definition 1. Given M is a set of matrices, define 


Pk{M, ||.||) := sup 




: Mi G M for 1 < i < fc 


where pk{M) is the largest possible norm of all products of k 
matrices chosen in the set M. The joint spectral radius p{M) 
is defined as 0251 


p{M) := lim pfc(M, ||.||)^. (7) 

k—foo 


independent of the induced norm used to define 

Fet us consider a set L of all possible adjacency matrices 
Ai and at each time point the adjacency matrix is randomly 
chosen from this set. L is surely bounded and may be finite 
or infinite. We define M as the set of system matrices corre¬ 
sponding to the adjacency matrices in L. Mi is a member of 
the set M and defined as Mi = [(1 — 5)1 + jdAi], Therefore, 
M is also bounded. If L is finite, M is also finite; if M is 
infinite, L is also infinite. 

Theorem 1. Consider a set L of all possible adjacency 
matrices of a dynamic network, with infection probability 
P, and recovery probability 5. If the joint spectral radius of 
set M of system matrices is less than one, the origin is an 
asymptotically stable equilibrium point and the epidemic dies 
out. 


Proof: Assume p(M) = I < 1. Definition [T] a value fco 
exists for each e > 0 with e < min(Z, 1 — 0 such that 

I — e < pk{M)'^ < I + e, yk > kg. 

If we raise all sides of the above inequality to the power k, 
we can conclude that 


lim p,,{M) = 0. 

k—foo 

Considering the formula of pk{M) in Definition [T] for any 
product of matrices Mi G M, we can write 


0 < 


n"< 


< pk{M). 


We have shown that the right hand side of the above inequality 
goes to zero when fc —oo. Therefore, for any product of MiS, 
we can write 


Therefore, the joint spectral radius of set M is the maximum 
possible norm of products of matrices in the set M when the 
number of products k goes to infinity. 


lim 

k—^oo 


riMi 


= 0 . 


Definition 2. Given M is a set of matrices, define 


Pk{M) := sup j P ( Mi j : Mi G M for 1 < / < fc > , 




where p denotes the spectral radius and 'pj.{M) is the largest 
possible spectral radius of all products of k matrices chosen 
in the set M. The generalized spectral radius 'p{M) is defined 
as lf26ll 


Considering the fact that ||A|| = 0 A = 0, if 


lim 

k—¥oo 


YlUMi 


= 0, we can write 


lim TT Mi = 0. 

k—foo ^ 


2=1 


If lim OiLi Mi = 0, for any initial infection probability 

k—foo 

vector Pq, we can write 


p{M) := lim p^{M)k _ (8) 

k—foo 

In 1271, the authors prove that, for a bounded set of matrices, 
the generalized spectral radius is equal to the joint spectral 
radius. 

Lemma 1. Four-member inequality. 1281 For a given arbi¬ 
trary set of matrices M and any fc > 1 

PkiM)^ < p{M) < p{M) < Pk{.M)^, 


lim 

k—^oo 


Rm. 


. 1=1 


Pq = lim 

k—¥oo 


Y[[{l-S)I + pAi 


.i=l 


Po=0, 


which shows that the origin is an asymptotically stable equi¬ 
librium point for any initial infection probability vector and 
any random sequence of adjacency matrices if the joint spectral 
radius is less than 1. In this case, the final infection probability 
vector is zero and the epidemic dies out. ■ 
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Theorem 2. Consider a set L of all possible adjacency 
matrices of a dynamic network with undirected graphs, set 
M of the system matrices corresponding to set L of the 
adjacency matrices, infection probability j3, and recovery 
probability 5. If the largest spectral radius of the matrices in 
set M is less than 1, the origin is an asymptotically stable 
equilibrium point and the epidemic dies out. 


Proof: If the network graph is undirected, its correspond¬ 
ing adjacency matrix is symmetric so that its corresponding 
system matrix is also symmetric. We know that for a given 
symmetric matrix Mi, we can calculate the induced 2 norm 
of Mi as 

\\Mi\\2 = ^Jp{Mi^Mi) = ^Jp{Mif = p{Mi). 


spectral radius. Hence, we can calculate the joint spectral 
radius as 


p{M) = lim = p((l - 5)1 + [iA). 

k—foo 

According to Theorem [T] the epidemic dies out if the joint 
spectral radius of the set of system matrices is less than 1. 
For a static network, the joint spectral radius is equal to 1 — 
5 + f3p{A). Therefore, the epidemic dies out if 


<5 < p{Ay 


( 11 ) 


The epidemic threshold for static networks in Eq. is 
identical to the one in its analytical results in m. 


If we use the induced 2 norm to calculate pi{M) and pi{M), 
we can write, respectively. 


Corollary 2. Consider a dynamic network with a fixed 
repetition pattern of T adjacency matrices in a set L = 
{Ai, A 2 , ■■■, At}, with infection probability (3 and recovery 


pi{M) = sup {WMiiy ■ S M} = sup{p{Mi) : Mi e M} probability 5. The epidemic dies out if 


and 


Pi{M) = sup{p{Mi) : Mi e M]. 

Therefore, we conclude that for a set of symmetric matrices. 

Pi (M) = Pi (M) = sup {p(Mi) : Mi e M} . (9) 

Moreover, we mentioned in Lemma [T] that the four-member 
inequality holds for any A: > 1. Therefore, we can write 

Pi(M) < p(M) < p(M) < pi(M). (10) 

Considering Eqs. (|9]) and (fTOl i. we can write 

p(M) = p(M) = sup{p{Mi) : Mi G M}. 

Therefore, the joint spectral radius of set M of symmetric 
matrices is equal to the largest spectral radius of matrices in 
the set. Based on Theorem [T] we conclude that if the largest 
spectral radius of system matrices of an undirected dynamic 
network is less than 1, the origin is an asymptotically stable 
equilibrium point and the epidemic dies out. ■ 

Corollary 1. Consider a static epidemic network with adja¬ 
cency matrix A, infection probability (3, and recovery proba¬ 
bility 5. The epidemic dies out if j < 

Proof: For a static network, M, the set of system matrices 
has only one element which is (1 — 5)1 -\- PA. In this case, 
p^(M), the largest possible spectral radius of all products of 
k matrices chosen in the set M, can be written as 


T 

p(n[(i-'))j+/3^i])<i- 

i=l 

Proof: Consider a dynamic network with a fixed repeti¬ 
tion pattern of T adjacency matrices and k = mT where m 
is a positive integer. For the case where k = mT, pj.(M) can 
be written as 


/ mT 


PkiM) =sup^p Mij I 
T 

n [(1 - s)i+pAi 

i^l 

V m 


= supjp 

/ T 


( 12 ) 


Because the set of system matrices is bounded, its joint 
spectral radius is equal to its generalized spectral radius. 
Hence, we can calculate the joint spectral radius as 


p(M) = lim PmriM) = P ( H + P^i] 

m—foo \ 

According to Theorem [T] if the joint spectral radius of 
the set of system matrices is less than 1, the epidemic 
dies out. Therefore, in this case, the epidemic dies out if 

p (^nLi [(1 “ ^ < 1 or equivalently 



Pk{M) = sup 



MiG M 


P{{1-5)I+PA)\ 


T 

p{Y[[{l-5)I + pAi])<l. 

2=1 


(13) 


M is a bounded set, and lIZTl proves that, for a bounded set The derived epidemic threshold for the periodic dynamic 

of matrices the joint spectral radius is equal to the generalized network is the same as the threshold in ||9|. 
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In the next corollary, we propose the condition under which 
epidemics die out in the case of dynamic networks with all 
the elements of set L corresponding to regular networks. In 
regular dynamic networks, although links between nodes are 
dynamic and a given node can change its neighbors at each 
index, all nodes have the same node degree, and they preserve 
their node degrees. 

Corollary 3. The epidemic in a dynamic regular network with 
undirected graphs and node degree of k dies out if ^ < k. 

Proof: We know that the spectral radius of a regular 
symmetric graph is equal to its node degree. Considering 
Mi = (1 — S)I + PAi, we can calculate the spectral radius 
of Mi for any regular adjacency matrix as follows; 


N 


^ El 


^q,n \ 


. 9=1 


= [l-d + {N- l)ppf Vn = 1,2,..., N, 

(15) 


where E denotes the expected value. Considering 


max, 

as 


\^q,n\ and (ffST l. we can calculate E < 
J = [l-,5 + (iV-l)^P]^ 


M 

|m 


Because the above equality holds for any product of MiS, 
E {pk{M)} can be written as 


e{MM)} = [i-s + {n-i)ppY 


— Consequently, the expected value of the joint spectral radius 

p{Mi) = 1 - ^ + Pp{Ai) = 1-6 +13 k. calculated as 


Therefore, all system matrices have the spectral radius of 
1 — 6 + pk, so the largest spectral radius of the system matrices 
is also equal to 1 — 5 + I3k. Based on Theorem |2] we can 
conclude that epidemics in dynamic networks with regular 
undirected graphs die out if the largest spectral radius of the 
system matrices is less than 1. Hence, the epidemic dies out 
ifl — 5 + /3fc<lor equivalently if 


E {p{M)} = E \ lim > = [1 — (5 + [N — l)f3P] ■ 

)-00 J 

We employ the Markov inequality and the expected value of 
the joint spectral radius to compute the upper bound for the 
probability that the joint spectral radius is more than 1. Using 
the Markov inequality, we can write 


Eq. (ITtIi remarkably is identical to the epidemic threshold for 
static regular networks 0. 

Hitherto, we have studied dynamic networks with deter¬ 
ministic and given adjacency matrices. However, for cases 
in which adjacency matrices are stochastic rather than de¬ 
terministic, the joint spectral radius of set of the system 
matrices is a random variable. Therefore, the condition under 
which epidemics die out is expressed in terms of statistical 
characteristics of the joint spectral radius. Dynamic Gilbert 
network is an example of such dynamic network with stochas¬ 
tic adjacency matrices. In a Gilbert network with parameter 
P, every link exists with probability P 1^ . In other words, 
the existence of a link is a Bernoulli random variable with 
parameter P. For a dynamic Gilbert network, at every time 
instant, these binary random variables are redrawn according 
to Bernoulli distribution. In the next corollary, we derive the 
upper bound for the probability of spreading epidemics for 
Gilbert dynamic networks in terms of the expected value of 
the joint spectral radius. 

Corollary 4. For a given dynamic Gilbert network with 
N nodes and probability P of the existence of links, the 
probability that an epidemic spreads is upper bounded by 
1 - 6 + {N - 1)I3P. 

Proof: We define M as M = 0^=1 where Mi 
denotes the system matrix corresponding to one realization of 
a Gilbert dynamic network’s adjacency matrix, fhq^n denotes 
the element in the row and column of matrix M. In 
the Appendix, we show that for all columns of M 


Prob{p{M) > 1) < E{p{M)}. (16) 

Substituting the expected value of the joint spectral radius in 
Eq. (O, we conclude that 

Prob{p{M)>l) <1-S+{N-1)PP. (17) 

According to Theorem [T] the epidemic dies out if the joint 
spectral radius is less than 1. Therefore, the probability of 
the epidemic spreading is equal to the probability that the 
joint spectral radius is greater than 1. Considering Eq. dll, 
we conclude that the probability of epidemic spread is upper 
bounded by 1 — i5 -f (TV — l)j3P. ■ 

However, when [1 — <5 -f {N — l)f3P] is greater than 1, 
this condition is not informative. Therefore, we consider the 
min{\, [1-5-1- {N — 1)/3P]} as the upper bound for the 
probability of spreading. 

IV. SIMULATION RESULTS 

In this section, we validate our theoretical results by sim¬ 
ulating an epidemic on synthetic and real dynamic networks. 
First, we simulate an epidemic on a dynamic Watts-Strogatz 
network and compare the derived epidemic threshold with 
the threshold proposed in Go). Then, the simulation result 
of the final fraction of infected nodes versus the joint spectral 
radius for a dynamic Barabasi-Albert network is presented. 
Further, we evaluate our analytical results for real networks 
by simulating an epidemic on the set of extracted graphs from 
the MIT Reality Mining data set 129)1 . Next, an epidemic on 
a dynamic regular network is simulated. Finally, we validate 
the derived upper bound for the probability of the epidemic 
spreading in a dynamic Gilbert network using the simulation 
results. 
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Fig. 3. Simulation results of epidemics on Watts-Strogatz and Barabasi-Albert dynamic networks (A) Final fraction of infected nodes for dynamic Watts- 
Strogatz networks with 1000 nodes and rewiring probability of 0.5. (B) Comparison between the joint spectral radius and the spectral radius of the product 
of system matrices for dynamic Watts-Strogatz networks. (C) Fraction of infected nodes over time for dynamic Barabasi-Albert network with 1000 nodes. 
(D) Final fraction of infected nodes for dynamic Barabasi-Albert network with 1000 nodes vs. the joint spectral radius. 


Fig. [3 illustrates the simulation results of epidemics on 
Watts-Strogatz and Barabasi-Albert dynamic networks. Both 
networks contain 1000 nodes, and the value of <5 = 0.2 
is kept constant while increasing /3. The final fraction of 
infected nodes is denoted by where I denotes the final 
number of infected nodes and N the total number of nodes. 
Fig. [3 (A) plots versus the joint spectral radius of system 
matrices and the spectral radius of the products of the system 
matrices for a dynamic Watts-Strogatz network with rewiring 
probability 0.5. In order to realize a dynamic Watts-Strogatz 
network, the adjacency matrix of the network at each index is 
chosen randomly from a set of four Watts-Strogatz graphs with 
average node degrees of 4, 8, 12, and 16 and, spectral radii, 
respectively, of 4.46242, 8.41081, 12.40911, and 16.38739. 
The joint spectral radius of the set of system matrices is equal 
to the spectral radius of the system matrix corresponding to the 
adjacency matrix with the largest spectral radius. We increase 
j3 from 0.00052 to 0.86652 in order to generate different 
epidemic strengths. The number of iterations for each case 
is 20. The observation is made that epidemics die out in all 
cases in which the value of the joint spectral radius is less than 
1. As soon as the value of the joint spectral radius increases 


beyond 1, the epidemic spreads, thus confirming the analytical 
results of Theorems [T] and |2] The curve of the final fraction 
of infected nodes versus the s^ctral radius of the system 
matrices product shows that p(ni=i [(1 “ 5)1 + is not 

an accurate epidemic threshold. Epidemics spread for some of 
its values that are less than 1. This simulation result contradicts 
analytical results of Theorem I in ifTOl that state that if the 
spectral radius of the product is less than 1, the epidemic dies 
out. 

Fig. [3 (B) compares the spectral radius of the product of 
system matrices of the dynamic Watts-Strogatz network with 
the joint spectral radius of its set of system matrices. The 
curve of Y = X allows us to determine the values of the joint 
spectral radius for which the spectral radius of the product 
is greater than the joint spectral radius or vice versa. Some 
points in Fig. [3 (B) have joint spectral radius greater than 1, 
and the spectral radius of the product less than 1, resulting in 
incorrect predictions of epidemics dying out for these cases if 
we choose the spectral radius of the product as the epidemic 
threshold. 

Fig. [3 (C) shows the fraction of infected nodes over time 
for a dynamic Barabasi-Albert network. In order to realize 
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Fig. 4. (A) Fraction of infected nodes for the MIT Reality Mining dynamic over time. (B) Final fraction of infected nodes for MIT Reality Mining dynamic 

vs. the joint spectral radius. (C) Final fraction of infected nodes for a dynamic regular network with 1000 nodes and node degree of 8. (D) Final fraction of 
infected nodes for a dynamic Gilbert network with 1000 nodes and node degree of 8. 


a dynamic Barabasi-Albert network, we have selected four 
Barabasi-Albert graphs with average node degree of 4, 8, 
12, and 16 with, spectral radii, respectively, of 12.66217, 
17.36462, 22.36887, and 27.91071, as the set of adjacency 
matrices. During simulation, the adjacency matrix of this 
dynamic network at each index is randomly chosen from this 
set of matrices with equal probability. The number of iterations 
for all cases is 20, and the initial fraction of infected nodes is 
0.2. Increasing the value of (3 leads to an increase in the joint 
spectral radius and, eventually, the final fraction of infected 
nodes. Also, the epidemics die out for cases in which the 
joint spectral radius is less than 1. Fig. [3 (D) depicts the 
final fraction of infected nodes for the dynamic Barabasi- 
Albert network where /3 increases from 0.00038 to 0.63481. 
Results show that epidemics spread for cases in which the 
joint spectral radius is greater than 1. 

Fig- El (A) shows the fraction of infected nodes over time for 
a dynamic network that has time-varying adjacency matrices 
extracted from the MIT Reality Mining data set ll29l . This data 
set contains the adjacency connectivity matrix of 94 persons, 
obtained by using mobile phones pre-installed with various 
specific software, including a logger of Bluetooth devices that 


was triggered when the distance between two mobile phones 
was approximately 5 m or less. Bluetooth scans were carried 
out every 5 min. This data set contains information collected 
from mobile phones from September 2004 to June 2005. We 
extract eight adjacency matrices for 8 consecutive hr, from 
8 : 00 a.m. to 4 : 00 p.m. on September 1, 2004. Specral radii 
of these matrices are, respectively, 6.30117, 5.41546, 9.44439, 
9.09696, 8.36535, 9.53451, 9.05251, and 7.41181. At each 
time index, the adjacency matrix was randomly chosen from 
the extracted matrices. For all the four cases, 5 = 0.2 and 
the initial fraction of infected nodes is 0.2. The number of 
iterations for each case is 50. Increasing /3 makes the joint 
spectral radius larger, and epidemic spreads for the values of 
the joint spectral radius larger than 1. 

Fig. |4] (B) depicts the final fraction of infected nodes for 
the dynamic network in Fig. |4] (A). We fix the value of <5 
at 0.2 and increase the value of /3 from 0.00089 to 0.83142. 
The number of iterations for each case is 50, therefore, the 
epidemic spreads when the joint spectral radius is greater than 
1. For cases in which the epidemic spreads, the final fraction of 
infected nodes increases with increasing value of joint spectral 
radius. 
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Fig. |4] (C) shows the hnal fraction of infected nodes for a 
dynamic regular network with 1000 nodes and node degree 
k = 8 versus the product of j and node degree. We 
simulate epidemics for various values of /3 in the interval 
[0.00106,0.99089], while hxing the value of S at 0.2. The 
number of iterations for each case is 100. Therefore, epidemics 
die out when j <P confirming the result of Corollary [3] 

TABLE I 

5 AND /3 VALUES USED IN SIMULATION OE THE EPIDEMIC IN THE 
DYNAMIC Gilbert network 


s 

/3 

1 - S + (N - 1)/3P 

0.95 

0.01 

0.09 

0.85 

0.01 

0.19 

0.74 

0.01 

0.3 

0.64 

0.01 

0.4 

0.54 

0.01 

0.5 

0.44 

0.01 

0.6 

0.34 

0.01 

0.7 

0.24 

0.01 

0.8 

0.14 

0.01 

0.9 

0.04 

0.01 

1 

0.7 

0.3 

1.5 

0.6 

0.4 

2 

0.5 

0.5 

2.5 

0.4 

0.6 

3 

0.3 

0.7 

3.5 

0.2 

0.8 

4 

0.1 

0.9 

4.5 

0.01 

0.99 

4.95 


Fig.|4](D) depicts the hnal fraction of infected nodes for a 
dynamic Gilbert network with 1000 nodes and a probability 
of connection P = 0.004 versus 1 — i5+ {N — 1)/3P, which is 
the upper bound for the probability of epidemic spread derived 
in Corollary |4| Table U shows the chosen values of 6 and /3 
in the simulation as well as values of 1 — <5 + {N — 1)I3P. 
The epidemic dies out up to the point at which the upper 
bound is less than 1. When this upper bound reaches 1, the 
epidemic begins to spread, conhrming the result of Corollary 
|4] Although having an upper bound greater than 1 for a 
probability is not informative, it can be used as a measure 
of epidemic strength. This just described situation in Fig. |4| 
(D), as increasing the value of the upper bound leads to an 
increase in the hnal fraction of infected nodes. 

V. CONCLUSION 

In this paper, we study the spread of SIS epidemics in 
dynamic networks. We propose an approach to derive the 
analytical epidemic threshold that can be applied to any 
dynamic network with an adjacency matrix randomly chosen 
from a set of matrices at each index. A linearized version 
of the nonlinear epidemic system is employed to derive the 
epidemic threshold. We show that an epidemic dies out if the 
origin is an asymptotically stable equilibrium point. We derive 
the epidemic threshold for dynamic networks using the joint 
spectral radius of system matrices. We calculate the simplihed 
version of the epidemic threshold for undirected dynamic 
networks based on the fact that the joint spectral radius of 
a set of symmetric matrices is equal to the largest spectral 


radius of the matrices in that set. Then, we derive the epidemic 
threshold for dynamic regular networks. For dynamic Gilbert 
networks, we compute the upper bound of the probability of 
an epidemic spreading in terms of the expected value of the 
joint spectral radius. 

Our analytical results show that epidemic thresholds of 
dynamic networks are determined by the dynamic topologies 
and epidemic strengths in the networks. In particular, the 
joint spectral radius of the set of system matrices determines 
whether or not the epidemic dies out. In other words, the 
joint spectral radius characterizes epidemic strengths and the 
level of connectivity between the nodes in dynamic networks 
over time. In the case of undirected networks, the joint 
spectral radius is dependent only on the adjacency matrix 
with the largest spectral radius, thereby implying that, for 
undirected networks, dynamics of epidemics are determined 
primarily by the adjacency matrix with the largest spectral 
radius. Variance in topology of dynamic networks impacts 
the spread of epidemics. In order to validate our analytical 
results for dynamic networks with various dynamic topologies, 
we simulate epidemics on Watts-Strogatz, Barabasi-Albert, 
Regular, MIT Reality Mining, and Gilbert dynamic networks. 
The dynamic Watts-Strogatz network is used to model a 
dynamic network with small-world properties, including a high 
clustering coefficient and small shortest path. On the other 
hand, the dynamic Gilbert network is employed to validate the 
derived epidemic threshold for dynamic networks with random 
structure and small clustering coefficients. Many networks in 
the real world, including e-mail networks, the world wide web, 
and biological networks, are considered scale-free networks 
that can be modeled by Barabasi-Albert networks EH-Ea. 
Our simulation result for dynamic Barabasi-Albert networks 
validates the derived epidemic threshold for such scale-free 
networks. Also, we verify our theoretical result with a real-life 
network. We simulate epidemics on the dynamic network with 
adjacency matrices extracted from the MIT Reality Mining 
data set 1291 . proving that the derived epidemic threshold holds 
for this real-life network as well. 

The model of dynamic switching networks in this paper 
does not take into account any temporal correlation between 
consecutive adjacency matrices. This assumption does not hold 
for dynamic networks with strong correlation between consec¬ 
utive adjacency matrices. For example, for the simulation of 
epidemics on the MIT Reality Mining data set, the adjacency 
matrix at each time step is chosen randomly and independently 
from other time steps. However, in reality, the sequence of 
adjacency matrices follows a temporal order, implying the 
existence of correlation between the adjacency matrices. The 
derived epidemic threshold corresponds to the most possible 
vulnerable sequence of matrices against epidemics, because in 
our framework, we assume that adjacency matrices are chosen 
completely at random and independent of each other. That is 
because the joint spectral radius of the set of system matrices 
is determined by the sequence of adjacency matrices that are 
the most vulnerable to epidemics. Thus, in the case of dynamic 
networks with temporal correlation, if the framework predicts 
the death of the epidemic, the epidemic will die out, even 
though this model ignores temporal correlations. However, the 
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prediction can be too conservative; That is, it may predict the 
spread of epidemic, when the epidemic actually dies out. 

Our analysis is based on the assumptions that the parameters 
and probabilities of the epidemics, the underlying structure, 
and topology of the transmission network are known. However, 
these assumptions will not hold in all scenarios. A robustness 
analysis with respect to model parameters can validate the 
application of the framework to such scenarios. We illustrate, 
all through this paper, some aspects of the robustness of 
our analysis by choosing adjacency matrices randomly and 
selecting widely ranging values of epidemic probabilities in 
the simulations. The framework developed in this paper may 
be extended to the derivation of epidemic thresholds for other 
types of epidemic models. It can be applied to cases in which 
the linearized version of the dynamical system corresponding 
to the epidemic model at the disease-free equilibrium point 
can be expressed in a matrix form similar to Eq. (01. 

Appendix 

Theorem 3. Consider M = Mi where Mi denotes the 
system matrix corresponding to one realization of a Gilbert 
dynamic network’s adjacency matrix. For the matrix M, the 
expected summation value of each column’s elements is 


= [l-d+{N- l)ppf Vn = 1, 2,..., N, 

as) 

where P denotes the probability of link existence and rh^ 
denotes the element in row and column of matrix M. 

Proof: We prove this theorem via induction. In the case of 
a Gilbert dynamic network, off-diagonal elements of adjacency 
matrix A are independent and identically distributed (iid) 
Bernoulli random variables with parameter P. The hrst step is 
to show that Eq. (fTsT i is correct when k = 1. Assume k = 1. 
In this case, M = (1 — S)I + (3A and E \^q,n 1} can 

be written as 



{ N ^ N-1 

= {l-6)+f3j2E{X,}, ( 19 ) 

q=l } i=l 

where the Xfs are iid random variables with parameter P. 
E{Xi\ = P. Therefore, we can rewrite Eq. (fT^ as 

= (l-<5) + (Af-l)/3-P. (20) 

Eq. (I 20 I 1 shows that Eq. (fTSl l is correct for k = 1 . The second 
step is to assume that Eq. (fTsT i is correct for k and prove it 
for k + 1. Assume M — Hti considering the assump¬ 
tion of correctness of Eq. m for k, = 

[1 — (5 -I- (TV — 1)/3P]^. Suppose R — Mk+iM. rq^n, the 
element in the the row and the column of R can 

be written in terms of the elements of M as 

N 

rq,n = (1 - S)mq^n + P ^ Xjfhj^n, (21) 


where the Xj's are iid Bernoulli random variables with pa¬ 
rameter P. Therefore, we can write as 


N 


N 


N 


\rq,n\ = (1 - (5) ^ fhq^n + I^Y 


N 


E 




( 22 ) 

where Xfs and Wj.n are independent. Hence, 


E{EYKn\} 


can be written as 


N 


N 


E S Vq.n\ f “ (f ^)E S T! f + 


>9=1 






N 


N 

E 


J=^,j¥=q 

However, E } 

1)^ Considering E 

[(1 — (5) -I- PP{N — 1)]^, we can rewrite (l23t as 


(23) 
(N - 


^ |e = [(1 - ^)+mN - 1)] E |e 


= [(1 - <5) +/3P(iV - 1)] 


fe+i 


(24) 


The result in Eq. (l24li for fc -f 1 is the last step in the proof of 
this theorem through induction. 
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